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The correlations among elements that break in random fuse network fracture are studied, for 
disorder strong enough to allow for volume damage before final failure. The growth of microfrac- 
tures is found to be uncorrelated above a lengthscale, that increases as the the final breakdown is 
approached. Since the fuse network strength decreases with sample size, asymptotically the process 
resembles more and more mean-field-like ("democratic fiber bundle") fracture. This is found from 
the microscopic dynamics of avalanches or microfractures, from a study of damage localization via 
entropy, and from the final damage profile. In particular, the last one is statistically constant, except 
exactly at the final crack zone (in contrast to recent results by Hansen et al., Phys. Rev. Lett. 90, 
045504 (2003)), in spite of the fact that the fracture surfaces are self-affme. 

PACS numbers: 62.20.Mk, 62.20.Fe, 05.40.-a, 81.40Np 



I. INTRODUCTION 

The scaling properties of fracture processes continue 
to attract interest from the statistical mechanics com- 
munity. Key quantities are the geometric properties of 
fracture surfaces and statistics of acoustic emission, or, 
in analogy to other systems, "crackling noise" . The point 
is that in failure of brittle materials the elastic energy of 
a sample is released in bursts. These "avalanches" often 
turn out to have scale-invariant statistics with respect 
to e i g. the probability distribution of the released energy 
Likewise, crack surfaces are often self-affme 
(with an empirical roughness exponent £) 0, Q, B| ■ The 
understanding of the origins of such critical-like statistics 
would perhaps be of interest to engineers ( "how to make 
tougher materials" ) but would also mean the solution of 
a very complicated many-particle system. 

In this respect, among the simplest models that are 
available are mean-field like fiber bundle models (FBM) 
UG3 and random fuse networks (RFN's) [Hill- The 
former describe democratic or global load sharing, and 
thus do not have anything close to the stress enhance- 
ments of real cracks (though one can introduce local load 
sharing to fiber bundles, and interpolate between these 
two limits as well). Such stress effects are to be found in 
a natural way in fuse networks, that simplify real elas- 
ticity by considering the electrostatic analogy. RFN's 
have two natural limits: weak disorder, when cracks are 
nucleated quickly and brittle failure takes place without 
much precursor activity, and strong disorder (without in- 
finitely strong elements ), w here damage develops before 
macroscopic failure UJ,UJ. 

The same signatures are found in the latter, RFN, case, 
that also characterize experimental systems: rough, self- 
affme cracks and microcracking that corresponds to the 
acoustic emission. The roughness exponent is in the prox- 
imity of Crfn ~ 0.7, in 2d, tantalizingly close to the 
minimum energy surface exponent, exactly 2/3. This re- 
sult holds also for e.g. 'weak' disord er [ll. [l5L [Tal and is 
close to what is seen in experimen ts Il7l(l 8l llfA l2(l | . The 
damage develops in avalanches @, Il0ll2l| . in analogy to 



democratic FBM's j2^|, so that one has for the probabil- 
ity distribution of the number of fuses (A) blown in one 
'event', for current-control, 

P(A)~A~ 5 / 2 . (1) 

The corresponding AE energy exponent is about [3 = 1.7 
[H |23| , as expected based on the exponent relation (3 — 
(5/2 + l)/2 I23. However, even for strong disorder finally 
stress enhancements come into play, and the sample fails 
catastrophically, with the elastic modulus (conductivity, 
in the RFN case) having a first-order drop. 

The purpose of this article is to investigate the devel- 
opment of the damage, between the MF-limit valid at 
the very initial stages of fracture and the final critical 
crack growth. Figure ^ shows an example of the transi- 
tion. The subplots depict the individual fuses that fail in 
subsequent parts of a stress-strain- (or current-voltage) 
history. Clearly, initially the damage is random (unless 
proven differently by more sophisticated analysis), and 
in the last panel it concentrates on the vicinity and at 
the final crack. 

In this respect, it is an important question how the 
pre-critical damage reflects the self-affine properties of 
the final fracture surface. Recently, Hansen and co- 
workers have attempted to relate its formation to a self- 
consistently developed damage profile that extends over 
all the sample 0, US ■ The scaling of the profile with 
the system size would then explain the roughness and 
its exponent. Clearly, this should also be visible in the 
dynamics of failure also prior to the end of the process. 
Another analogy is given by dynamics in dipolar ran- 
dom field magnets, which can account for the symmetry 
breaking (as signaled by the formation of the final crack) 
due to shielding in the direction of the external voltage 
and for stress enhancements that drive cracks mostly per- 
pendicular to it |2q. 

We study these aspects by concentrating on two kinds 
of quantities: those that characterize the spatial distri- 
bution of damage in samples, and those that analyze the 
temporal correlations in individual failure events (as e.g. 
during an avalanche, or series of fuse failures due to the 
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FIG. 1: Snapshots of subsequent damage patters in a failure 
of a RFN (each sub-plot having the same number of failed 
fuses, separately). R — 1, "strong disorder", L = 60. 



increase of a control parameter). Section II considers the 
former, and uses as the main tool entropy, comparing the 
damage integrated over windows of time and/or space to 
that in, spatially, completely random damage formation. 
We also study the damage profiles of completely failed 
samples. From both kinds of analysis emerges a picture 
of crack development, that is mean-field -like beyond a 
finite interaction range, and until the final breakdown is 
induced by rare event statistics Jl^. This in particu- 
lar includes the fact that except in the "fracture process 
zone" , i.e. in the vicinity of the final crack, the damage 
is statistically homogeneous. Thus in this particular case 
of RFN's the theory proposed by Hansen et al. seems un- 
likely to be the explanation for the self-affine geometry 
of cracks. 

In Section III, the internal dynamics of avalanches is 
considered. We look at the probability distributions and 
average values of "jumps" (relative changes in the po- 
sition of subsequent failures). It transpires that there 
is a smooth development, in which the these quantities 
exhibit a cross-over from the FBM/MF-like lack of spa- 
tial inhomogeneity towards localized crack growth within 
a scale £. This resembles some observations by Curtin 



about critical damage clusters in a more elaborate fiber 
bundle-type model: the dynamics is based on democratic 
load sharing in spite of the presence of stress enhance- 
ments [27|]. It also pertains to the question of the exis- 
tence of "representative volume elements" [2^ or coarse- 
graining in fuse networks [23] , related to the general ques- 
tion of how to account for microscopic dynamics and phe- 
nomena with coarse-grained variables and equations. Be- 
yond any such correlation length £ as may exist within 
avalanches the network looks homogeneous, if in addition 
the damage density is statistically homogeneous. Finally, 
we finish the paper in Sec. IV with a summary and some 
open prospects. 



II. DISTRIBUTION OF DAMAGE 

The RFN's, as electrical analogues of (quasi-)brittle 
media consist of fuses with a linear voltage-current re- 
lationship until a breakdown current ib- A stress-strain 
test can be done by using adiabatic fracture iterations: 
the current balance is solved, and at each round the most 
strained fuse is chosen according to min(ij/j CJ ), where 
ij is the local current and j c j the local threshold) . Cur- 
rents and voltages are solved by the conjugate-gradient 
method. In the following we use for P(i c ) a flat distribu- 
tion P(i c ) = [1 — R, 1 + R], with the disorder parameter 
R chosen as unity. The simulations are done in 2d, in 
the (10) lattice orientation, with periodic boundary con- 
ditions in the transverse direction (y). Square systems 
upto 100 2 have been studied; notice that the damage is 
in practice volume-like, and thus thousands of iterations 
are needed per a single system for L ~ 100. 

Studies of the break-down current lb as a function 
affirm the expected outcome of a logarithmic scaling 
[r2l[l^ |. resulting from extremal statistics (lb ~ L/hxL). 

T 2 

This implies in the mean-field limit that rib ~ logL , where 
rib is the average number of broken fuses in a system. 

To analyze the spatial of distribution of damage it is 
useful first to take note of the fact that in the latter stages 
the system behaves anisotropically: just before the for- 
mation of a critical crack the spatial density of the bro- 
ken fuses should be a stochastic variable, with a constant 
mean in the transverse direction to the external voltage. 
However, along the voltage direction differences may en- 
sue. To study such trends in the damage mechanics and 
the localization we consider the entropy of the damage 
averaged over y in each sample (this is in analogy to the 
procedure used with AE experiments of Guarino et al. 
[2j). The network is divided along the current flow direc- 
tion into sections, and the entropy S defined as 



S 



(2) 



where qi is the fraction of burned fuses in section i. S is 
normalized by S e , the entropy of a random, on the aver- 
age homogeneous distribution of failures (of equal total 



damage). Thus the extreme limits are zero and unity, 
corresponding to completely localized damage and com- 
plete random one, respectively. The final crack extends 
between y — and y — L, and a sensible choice is to use 
for the section width Sx a value larger than the typical 
interface width w, 



w=({h y -hf) 1 '\ 



(3) 



where h y denotes the crack location and h its mean po- 
sition in the x-direction. Since w ~ ZA, with f < 1, it is 
clear that using a constant number of sections will with 
increasing L localize the fracture zone either entirely in- 
side one, or between two neighboring ones. For better 
statistics it is preferable to have Sx >> 1, though the 
interpretation is perhaps more difficult than for the ex- 
treme value 5x = 1, say. The width of the sections used 
in computing the entropies was chosen to be 6x = L/W. 
In this discrete form, the entropy reads 
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FIG. 2: S vs. l/(LlogL) in networks with strong disorder 
(7? = 1). The straight line is a linear least squares fit that 
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(4) 



where k = L/Sx is the number of sections, rij the num- 
ber of burnt fuses in the «'th section, and N — ^Zj=i n i 
is again the total. Note that the absolute value of S is 
dependent on the choice for Sx. S can now be used to 
consider different parts of the stress-strain curve, sepa- 
rately, or the final damage pattern. 

Fig. El shows the total entropy versus system size. 
The best kind of linearity with regards the data is ob- 
tained with a scaling variable 1/(L log L) . This can be 
considered with the following Ansatz. Assume, that 
the fractures are distributed otherwise randomly (n^), 
except the one containing the final crack, which has 
n;,, fc + An^fc. Take An;, <C (nb), which implies approxi- 
mately S w 1 — \nk(Arib.k/ n b)- Noticing the logarithmic 
scaling of damage, it follows that S ~ — l/(£logL), if 
and only if Arib,k scales as An^k ~ L/(logL) 2 . We have 
not checked explicitly that this holds; note that the frac- 
ture surface, being self-affine, is then supposed to contain 
L a fuses, with 1 < a « 2, but the fracture process zone 
itself, contains other damage (broken fuses) contributing 
to Anb (see Fig. 1 again, and the last panel in particu- 
lar). In any case, it is obvious that the entropy increases 
with system size, indicating more and more completely 
random damage. The limiting value of S is slightly below 
unity; it is hard to say whether this difference is due to 
the choice used for computing S or a real one. 

In figure |3] an example is shown of how the damage 
actually localizes when the fracture process is divided 
into sequential slices. It is clear that initially most of 
the fuses break randomly, and only in the last one strong 
localization takes place. Similar sampling can be done 
also with the external voltage or current as control pa- 
rameters, the difference between these two being that the 
fracture is more abrupt in the latter case. 

This lack of the localization of damage is reflected in 
Fig. E] The damage density (p)(x) has been averaged 
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FIG. 3: Entropy S versus time interval An;. L — 100, R = 1. 
Average over 20 realizations. 



in the y-direction as a function of the normalized coor- 
dinate transform \— [x — x c + L)/2L, where the point 
x c is chosen as the one with the maximum damage, and 
is located in practice at the final fracture line, x c ~ h. 
After this shift, the average density is computed taking 
care that it is normalized correctly since the number of 
samples contributing for each \ varies with the final crack 
location - x c is a random variable. We also have added 
the average fracture line width w(L = 100) as a compar- 
ison (from ref. 16]). It can be seen that outside of the 
immediate vicinity of the fracture process zone the dam- 
age is constant. Notice the error bars of the data points, 
and that the data points located far away from the crack 
line suffer from the presence of less data points as seen 
from the error bars. It would be interesting to analyze 
in detail the functional shape of (p)(x) i n the proximity 
of the crackline, x = 0-5- The implication of the results 
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FIG. 4: Averaged damage density (p)(x)- L = 100, R=l. 



between consecutive fractures, 

Ax = \x i+ i - Xi\ 

&y = \vi+i - ml 



(•5) 
(6) 



where Xi and t/j are the x- and y-coordinates of the ith 
fracture. Another one choice is given by the average 
distances between consecutive fractures belonging to the 
same avalanche (ie. induced by a single increment of the 
control parameter) 



A-l 



avalanche 



1=1 

1 A_1 

Ay ava i anc fi e — — — ^ ^ |yz+i 



1 - Xi 



(7) 
(8) 



is that the density can be written as a sum of a constant 
(L-dependent) background, and a term that has to decay 
(perhaps exponentially) within a finite lcngthscalc from 
the crack. This decay length in turn may depend on L. 

Such an observation is in contradiction to the proposed 
"self-consistent" quadratic functional form, by Hansen et 

al. Q. This would imply {p)( X ) = p f - A ( L(2 ; X x ~ 1) ) , 
which is clearly not the case. In the light of the picture 
discussed below about the internal dynamics of microc- 
racks or avalanches, the interpretation is that the final 
crack is formed here similarly to weak disorder in a "crit- 
ical" manner. That is, once a damage density sufficient 
for "nucleation" is established the largest crack becomes 
unstable. Prior to that the correlations in the damage ac- 
cumulated can for all purposes be neglected. This would 
in turn to imply that the origin of the self-afnne crack 
roughness in fuse networks is not dependent on whether 
there is "strong" or "weak" disorder, as long as there 
are no infinitely strong fuses, or as long as the process 
does not resemble e.g. percolation due to the complete 
domination of zero-strength fuses. 



III. AVALANCHES 

Next we consider the correlations in the dynamics of 
individual fuse failures. Recall that the MF- limit states 
that consecutive ones should not be spatially correlated 
in any fashion; the opposite limit is given by the growth 
of a linear crack in which it is always the one adjacent to 
the crack tip to fail next. In the case that the growing 
crack is "rough" one expects that the subsequent failure 
takes place inside a fracture process zone, analogously to 
normal fracture mechanics, one of the follow-up questions 
being how the size and the shape of this zone vary with 
system size and as the crack grows |l5j . The simplest 
quantities to compute, to examine localization and spa- 
tial correlations between fractures, are the Id distances 



The MF theory predicts P(Ax = k) = 2(L - k)/L 2 and 
P(Ay — k) = 2/L, if the boundary conditions used here 
are taken into account, and A is again the avalanche size 
measured in the number of fuses broken during it. 

Fig. [3] depicts, as a comparison for the mean-field re- 
sults, the average distances in the x- and y-dimensions 
between consecutive broken fuses belonging to the same 
avalanche. (Ax ava i anche ) and (Ay avalanche ), are shown, 
respectively, as a function of the system size. Both are 
linear like in the MF theory, but with a smaller slope 
with L. This means, that the damage created by a typical 
avalanche (microcrack creation, crack advance etc.) is lo- 
calized compared to the MF-prediction, but nevertheless 
the localization does not get stronger with L. One should 
note that the damage as such is almost volume-like. The 
result is thus not surprising in the sense that a reduction 
of the slope (sublinear behavior, say, (Ax ava ianche) ~ L a , 
with a < 1) would imply concomitant faster average 
crack growth, which would be in contradiction with the 
damage scaling. 

To understand in detail the dynamics of microcracks 
is a difficult task. This is since the growth dynamics 
is not local: the burned fuses do not have to form con- 
nected clusters by any remotely easy criterion. It is easy 
to comprehend that the driving force for the localization 
is standard stress-enhancement, but as is true for RFN's 
crack shielding and arrest (due to strong fuses, in the 
early stages of fracture) play a role. One may set aside 
for the sake of discussion the separation of the events 
into avalanches, and just consider the distances between 
consecutive burned fuses. Figure demonstrates the dif- 
ference between two probability distributions P(Ax), av- 
eraged over the first 1/8 of the typical failure process and 
the last, respectively, for a fixed L. As one could expect, 
there is a peak in the distribution (this holds for both x- 
and y-directions separately), which is greatly suppressed 
in the first part closest to the MF- limit. 

Thus one may conclude, that there is a continuous 
cross-over from purely MF-like behavior to a complicated 
non-local growth dynamics. This is also exhibited by 
such distributions P(Ax), P(Ay). The analysis of the 
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FIG. 5: Average one-dimensional distances (Ax avalanche) and 
(Ayavaianche) between consecutive broken fuses belonging to 
the same avalanche vs. linear dimension of system L. R = 1. 
The solid lines are linear least squares fits with slopes 0.018 
for (Axavaianche) and 0.055 for (Ay avalanche) ■ The dashed 
lines, linear with slopes 0.33 for (Ax avalanche) and 0.25 for 
(Ay a vaianche) , correspond to the distances predicted by mean- 
field theory. Statistical errors are smaller than or equal to the 
size of the data points. 
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FIG. 6: Distribution of average distances in the a;-dimension 
between consecutive broken fuses computed over the 1 st and 
8 th 8 th of the fractures of 20 realizations. L = 100, R = 1. 
The solid lines correspond to mean-field predictions. 
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FIG. 7: The scaling of £ H (see Eq. ©) with increasing damage 
for L = 100. The total number of broken bonds has been 
divided into ten consecutive windows, and in each of these 
(d y ) has been computed, and £ y using Eq. ©. For r two 
guesses (1/2, 2/3) are used, note that < r < 1. 



correlation length), £ x or £ B . Second, one should recall 
the scaling of the strength with L: catastrophic crack 
growth takes place earlier and earlier with respect to the 
intensive variable, current. This means that the RFN's 
resemble, in the thermodynamic limit, more and more 
the mean-field-case in their fracture properties in spite 
of the stress- (or more exactly current-) enhancements 
that the model contains. 

Again, our data does not allow us to conclude firmly 
how such correlation lengths behave for Ax or Ay small - 
how the associated distributions P would scale for small 
arguments that is. One may however simply use an 
Ansatz that P ~ Ay~ T upto £ y , say and MF-like for 
larger Ay This defines the correlation length £ y for 
a given damage density p. Using now the distribution P 
allows to compute (Ay) and relate it to £ y , valid for such 
deviations from the uncorrelated fracture process, that 
£ v <g (L/4). The result is in analogy to Delaplace et al. 
that 



detailed shape of the small- argument part of the proba- 
bility distributions would be an interesting challenge. To 
first order, the result is a convolution of a microcrack 
size distribution and the corresponding stress enhance- 
ment factor, such that the distribution P evolves accord- 
ing to the growth indicated by the involved probability 
distributions. Given the simple forms of say P(Ax) for 
small arguments there might be some hope for developing 
analytical arguments. 

When considered as a function of L it becomes imme- 
diately apparent why the avalanche statistics resembles 
the MF-case so much. This is due to two separate factors: 
first, the growth is clearly in the sample case of Fig. 
(or Fig. 1 again) local over a certain lengthscale (damage 



(dy) = 



1-t {2-T)L 2 +8ri 



4(2 -r) (1-t)L + 2t£ v 



(9) 



In the opposite limit, £ y ~ 1 the correlations are badly 
defined since the model is discrete. Figure [3 shows an 
example of the ensuing scaling with different guesses for 
t, for L — 100. The main observation is an exponential 
(perhaps) increase of t; y with damage. Again, note that 
with still larger system sizes the total damage is dimin- 
ished, which in turn implies that the maximal correla- 
tions in the damage accumulation become weaker. Please 
observe that we have not studied in detail the other possi- 
bility, £ a since the main interest lies with the correlations 
in the crack-growth direction. 
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IV. SUMMARY 

In this article we have studied the distributions and 
development of damage in random fuse networks, with 
"strong" disorder. Our aim has been to understand pos- 
sible deviations from mean-field theory, and the associ- 
ated correlations. This is of relevance both as regards the 
statistical mechanics of fracture in general, and in partic- 
ular also the growth and formation of self-afhne cracks. 
In other words, we have concentrated on the "approach 
to the critical point" if the failure transition is considered 
as an analogue of ordinary phase transitions. 

An analysis of the localization of damage both during 
the fracture process and a posteriori reveals that in the 
case studied the correlations are very weak, are formed 
mostly in the last catastrophic phase of network failure, 
after the maximum current I ma x, and do not have any 
global correlations. The localization is centered in and 
around the final crack surface, or what may be called 
as the total volume encompassed by a "fracture process 
zone". We would like to note that this is in contrast to 
the recent theory of self-organized damage percolation, 
of Hansen et al. (|24j) devised to explain the formation 
of self-affine cracks in fuse networks, and in related ex- 
periments. In particular it should be stressed that there 
is no evidence of a global, non-trivial damage profile as 
contained in that proposition. Recent numerical results 
of Nukala et al. |30| with much better statistics than 
what is the case here or with other, earlier authors also 
seem to imply the same. 

Since also in this particular case much of the dam- 
age incorporated in the final crack is due to the "last" 
avalanche it seems then logical that the fracture surface 
geometry is formed similarly to RFN's with weak dis- 
order, for which there is no quasi-volume like damage, 



and one often has just the propagation, and formation 
of the final crack. Nevertheless of this complication, the 
measured roughness exponents from such simulations are 
close to that in the case here at hand, but note also that 
such exponents are notoriously hard to measure numeri- 
cally in finite sized samples: more extensive work in this 
respect would certainly be desirable. 

The internal correlations of the avalanches become 
more and more important as damage grows, but in line 
with the fact that the statistics is close to the mean- 
field case the growth is never very far from the MF: for 
all phases studied there are remote, broken fuses instead 
of ones localized close to the last growth event. There 
is an associated lengthscale that can be roughly defined 
based on the x- and y-dependent results, but of course 
one could go further and look at the radial probability 
distribution P(r), with r = (x i+ i,y l+1 ) - {x^yCj (for 
which one would presumably need still much larger sys- 
tems to get decent averaging) . One central lesson is that 
localization will diminish with system size due to the nor- 
mal volume effect of strength, decreasing with L. In this 
respect, fuse networks are not unique, and other simu- 
lation models of brittle fracture should exhibit the same 
behavior. To conclude, even in our case with quite strong 
disorder the failure process consists of weakly correlated 
damage growth and a final catastrophic crack propaga- 
tion phase, that induces a first-order drop in the elastic 
modulus. 
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